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MAXIMAL LP ANALYSIS OF FINITE ELEMENT SOLUTIONS 
FOR PARABOLIC EQUATIONS WITH 
NONSMOOTH COEFFICIENTS IN CONVEX POLYHEDRA 

BUYANG LI AND WEIWEI SUN 


Abstract. The paper is concerned with Galerkin finite element solutions of 
parabolic equations in a convex polygon or polyhedron with a diffusion co¬ 
efficient in W 1 ' N+a for some a > 0, where N denotes the dimension of the 
domain. We prove the analyticity of the semigroup generated by the discrete 
elliptic operator, the discrete maximal L p regularity and the optimal L p error 
estimate of the finite element solution for the parabolic equation. 


1. Introduction 

Let 17 be a bounded domain in (with TV = 2 or TV = 3), and let Sh be a 
finite element subspace of Hq (17) consisting of continuous piecewise polynomials of 
degree r > 1 subject to certain quasi-uniform triangulation of the domain 17. We 
consider the parabolic equation 

( dtu — V • (aVu) = / in 17 x (0, oo), 

(1.1) < u = 0 on <917 x (0, oo), 

[ u(-, 0) = u° in 17, 

and its finite element approximation 

n 9 n, / (■ d t u h ,v h ) + (aVu h , Vv h ) = (f,v h ), V v h e S h , 

{ ’ Wo) = <, 

where / is a given function, and a = ( dij(x))N X N is an TV x N symmetric matrix 
which satisfies the ellipticity condition 

N 

(1.3) A _1 |£| 2 < ^ aij(x)€i£j < A|£| 2 , for x G fl, 

i,j =1 

for some positive constant A. 

If we define the elliptic operator A : Ro(f7) —>• i7 _1 (f7) and its finite element 
approximation Ah : Sh —t Sh by 

(1.4) (Aw,v) := (aVw,Vv), V w, v G Hq (17), 

(1-5) (A h w h ,v h ) := (dVw h ,Vv h ), V w h , v h £ S h , 
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then the solutions of EB and can be expressed by 


( 1 . 6 ) 

(1.7) 


u(t) = E(t)u° + f E(t — s)/(s)ds, 

J o 

u h (t) = E h (t)u° h + [ E h (t - s)/(s)ds, 

J o 


where {E(t) = e~ tA }t>o and {Eh(t) = e~ tAh }t >o denote the semigroups generated 
by the operators —A and — Ah, respectively. By the theory of parabolic equa¬ 
tions and [32], it is well known that {E(t)}t>o is an analytic semigroup on Co(f2) 
satisfying 


(1.8) \\E(t)v\\ L °°+t\\d t E(t)v\\ L °° <C\\v\\ L ™, Vv e Co(n), Vt > 0, 


which is equivalent to the resolvent estimate 

||(A + A)- 1 i;|| L o 0 KCX^Mlo., VveCo(n), VAeS 8+l/2l 

where Sg + w /2 := {z £ C : |arg(z)| < 9 + 7r/2}. The counterparts of these two 
inequalities above for the discrete finite element operator Ah are the analyticity of 
the semigroup {Eh(t)}t> o on L°° D Sh- 

(1.9) \\E h (t)v h \\ L °° +t\\d t E h {t)vh\\L°° < C\\v h \\ L °°, Vv h G S h , Vt > 0, 


and the resolvent estimate 


II (A + Ah) 1 h/i||l°° < C\ 1 \\vh ||l°° j V Vh & Sh, V As ^ v +n/2- 

The estimates of the discrete semigroup have attracted much attention in the past 
several decades. With these estimates, one may reach more precise analyses of finite 
element solutions, such as maximum-norm analysis of FEMs ED S3 m EH] , error 
estimates of fully discrete FEMs pH 132 m2 and the discrete maximal L p regularity 
for parabolic finite element equations pj E3 EH OH HU- 

The proof of (11.911 dates back to Schatz et. al. [35], who proved (11.91) with a 
logarithmic factor for the heat equation in a two-dimensional smooth convex domain 
with the linear finite element method. The logarithmic factor was removed in the 
case r > 4 for TV = 1, 2,3 in [32], and the analysis was further extended to the case 
1 < N < 5 in [J. Later, a unified approach was presented in [32] by Schatz et. al., 
where they proved m with the Neumann boundary condition for all r > 1 and 
N > 1. The result was extended to the Dirichlet boundary condition in m for the 
linear finite element method. Some other maximum-norm error estimates can be 
found in Ei e nm Els m na , and the resolvent estimates can be found in mm- 
A related topic is the discrete maximal L p regularity (when u° = 0 and 1 < 
p,q < oo) 

(1-10) \\dtUh\\LP((0,T)-,LQ) + \\AhUh\\LP((0,T)-,LQ) < II / \\lp((0,T)-,Li) ; 

which resembles the maximal L p regularity of the continuous parabolic problem 
and was proved by Geissert mug. A straightforward application of (11.101) is the 
L p - norm error estimate 


(1.11) || PhU — Uh\\LP{(0,T)-,L*) < C Pt q(\\PhU° — Uh\\Li + II PhU — Rhu\\LP((0,T)-,Li)), 

where Rh is the Ritz projection associated with the operator A and Ph is the L 2 
projection onto the finite element space. 
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All these estimates were established under the assumption that the coefficients 
aij and the domain fl are smooth enough so that the parabolic Green’s function 
satisfies 
( 1 . 12 ) 

\dfd^G(t,x,y)\ < C(t 1 / 2 + |a - |/|) _(JV+27+lfll) e _]£ ^, V 0 < 7 < 2, 0 < |/ 3 | < 2 . 


Although the condition on the coefficients was relaxed to £ C' 2+ “(U) in M, this 
assumption is still too strong for many physical applications. One of the examples 
is an incompressible miscible flow in porous media EE EH, where the diffusion- 
dispersion tensor [ajj]fj =1 is only a Lipschitz continuous function of the velocity 
held. In a recent work [2?] , the first author proved 0 in a smooth domain under 
the assumption aij £ kF 1,00 (f2), together with the estimate (when u° = 0 and 
1 < p, q < 00) 

(1-13) II^/i||l p ((0,T);VU 1 ’ , z) < C’p t q\\f\\ L p^ 0tT y iW -i, q ' ) , 

which were then applied to the incompressible miscible flow in porous media m- 
Moreover, the problem in a polygon or a polyhedron is of high interest in practical 
cases, while the inequality (11.121) does not hold in arbitrary convex polygons or 
polyhedra, and all the analyses of (11.1011 - (11.1311 are limited to smooth domains so 
far. For the problem in two-dimensional polygons with constant coefficients, the 
inequality 0 with an extra logarithmic factor was proved in mmim by using 
the following estimate of the discrete Green’s function T^: 


\Th(t,x,Xo)\dx < C | ln/i|. 


The corresponding results in three-dimensional polyhedra are unknown. More 
interested is whether these stability estimates hold with the natural regularity 
£ TU 1,P (U) for some 1 < p < 00, since such estimates are important for the 
extension of the analysis to a general nonlinear model. 

This paper focuses on (II,911 - (II.1011 and (11.131) in a convex polygon or polyhe¬ 
dron with a weaker regularity of the diffusion coefficient. Instead of estimating T^ 
directly, we present a more precise estimate for the error function F := T^ — T 
(see Lemma 2.2) with which the logarithmic factor can be removed (this idea was 
used in [25]), where T is a regularized Green’s function. To compensate the lack of 
pointwise estimate of the second-order derivatives of the Green’s function, we use 
local W 1 ' 00 estimate and local energy estimates of the second-order derivatives (see 
Lemma 4.1). Our main result is the following theorem. 


Theorem 1.1. Assume that £ W 1,N+a {0) for some a > 0, satisfying the 
condition 0, and assume that Q is either a convex polygon in R 2 or a convex 
polyhedron in R 3 . Then 

(1) the semigroup estimate (11.91) holds, 

(2) the solution of (|1.2D satisfies Ijl.lOD when f £ L p ((0,T);L 9 ) and u° = 0, 

(3) the solution of (11.21) satisfies (11.131) when f £ L p ((0,T); W~ 1,q ) and u° = 0. 


Under the assumptions in Theorem ll.il and assuming that the solution of 0 
satisfies u £ C(Q x [0,T]), (11.111) follows immediately from (11.101) . 

The rest of this paper is organized as follows. In section 2, we introduce some 
notations and present a key lemma based on which our main theorem can be proved. 
In section 3, we present superapproximation results for smoothly truncated finite 
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element functions and present several estimates for the parabolic Green’s functions 
under the assumed regularity of the coefficients and the domain. Based on these 
estimates, we prove our key lemma in section 4. 


2. Notations, assumptions and sketch of the proof 


2.1. Notations. For any nonnegative integer k and 1 < p < oo, we let W k,p (Q) 
be the conventional Sobolev space of functions defined in f2, and let Wg 1,p (f2) be 
the subspace of kF 1,p (fi) consisting of functions whose traces vanish on dfl. As 
conventions, we denote the dual space of W 0 1,p (fi) by W~ 1,p (f l) for 1 < p < oo, and 
denote H k (fl) := W k ’ 2 (Q) and L p ( fl) := VF°’ p (fi) for any integer k and 1 < p < oo. 

Let Qt := fi x (0,T). For any Banach space X and a given T > 0, we let 
L p ((0, T); X) be the Bochner spaces equipped with the norm 


||/||i*>((0,T);X) 


( j o ll/Wllx^ > 

ess sup ||/(i)||jr. 
te(o,T) 


1 < p < 00, 

P= oo, 


To simplify notations, in the following sections, we write L p , H k and W k,p as the 
abbreviations of L p (fl), H k (Q) and W k ’ p (VL), respectively, and denote by (•, •) the 
inner product in L 2 (fl). For any subdomain Q C Qt , we define 

Q l := {a; G : (x,t) G Q}, 


il/llz,°°' 2 (Q) :=ess sup ||/(-,*)||z,= ( q*), 
te(o,T) 

\\fhp(Q) ■= (^JJ Q \f(x,t)\ p dxdt^J , VI < p < oo, 

and denote w(t) = w(-,t) for any function w defined on Qt- 

We assume that f l is partitioned into quasi-uniform triangular elements T k . I = 
l,--- ,L , with h = max;{diamr/'}, and let Sh be a finite element subspace of 
consisting of continuous piecewise polynomials of degree r > 1 subject to 
the triangulation. Let a{x) = (aij(x)) NxN be the coefficient matrix and define the 
operators 


by 


o T ±1 i Ah : Sh —t Sh, 


A: HiH~ l 

Rh-.Hl^Sh, Ph-L 


S h , 


(A<j>,v) = (aV0, Vu) 
(A h <t> h ,v) = (aX(f> h , Vw) 
(• A h R h w,v ) = (Aw,v) 
(P h <P,v) = (<j),v) 


for all </),v G Hi , 
for all 4>h G Sh, v G S h , 
for all w G Hi and v G Sh, 
for all </)Gi 2 and v G Sh- 


Clearly, Rh is the Ritz projection operator associated to the elliptic operator A 
and Ph is the L 2 projection operator onto the finite element space. The following 
estimates are useful in this paper. 
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Lemma 2.1. If Q is a bounded convex domain and aij £ !L 1 ’ Ar+Q (f2), N >2, then 
we have 

(2.1) IMIff 2 < C||v- (aVu;)|| L 2 , \/w £ tf 0 \ 

(2.2) ||Vui||loo < Cp||V • (oVu;)||lp, for any given p > N , \/w£Hq, 

and the solution of CUD with u° = 0 satisfies 

(2.3) H<9iU||LP((0,T);L<!) + ||^4u|| iP (( 0)r ) ;i ,q < C P: q || /1| lp((0,T)-,LQ) > 

(2.4) \\dtu\\ L p(( 0y T);W- 1 ’i) + IM|lp(( 0 ,T);IU 1 ..j) < Cp ig ||/||LP(( 0i T);IU-U9), 

for all 1 < p, q < oo. 

In the Lemma above, (EH) is the standard U 2 -regularity estimate in convex 
domains and (|2.2D is a simple consequence of the Green’s function estimates given 
in Theorem 3.3-3.4 of [18], and (I2.3I) - (I2.4I) are consequences of the maximal L p 
regularity (see Appendix for details). 

2.2. Properties of the finite element space and Green’s functions. For any 

subdomain D C Cl, we denote by Sh(D) the space of functions restricted to the 
domain D, and denote by S®(D) the subspace of Sh(D) consisting of functions 
which equal zero outside D. For any given subset D C Cl, we denote Bd{D) = {x £ 
fl : dist(a;,H) < d} for d > 0. Then there exist positive constants K and k such 
that the triangulation and the corresponding finite element space Sh possess the 
following properties ( K and n are independent of the subset D and h). 

(PO) Quasi-uniformity: 

For all triangles (or tetrahedron) r[ l in the partition, the diameter hi of r/ 1 and 
the radius pi of its inscribed ball satisfy 

K~ 1 h < pi < hi < Kh. 

(PI) Inverse inequality: 

If D is a union of elements in the partition, then 

\\Xh\\w‘.p(D) < Kh N ^ p ' > \\Xh\\w k ’i{D), V Xh £ Sh, 

for 0 < k < l < 1 and 1 < q < p < oo. 

(P2) Local approximation and superapproximation: 

(1) There exists a linear operator A, : i?g(f2) —> Sh such that if d > nh, then 

k 

\\v-I h v\\ L 2 (D) <K^h k d~ l \\v\\ H k-i {Bd{D)) , V v £ H k n H^, 1 < k < 2. 

;=o 

Moreover, if supp(u) C D , then IhV £ Sfl(Bd(D)). For example, the Clement 
interpolation operator defined in [5] has these properties. Also, the Lagrange inter¬ 
polation operator 11^ satisfies 

11^ - Rhv\\ L 2 {D) + h\\V(v - U h v)\\ L 2 {D) < A7i 2 ||V 2 u|| L 2 (Bd(D)) , V v £ H 2 n H], 

(2) If d > nh , uj = 0 outside B 2 d(D) and Id^wl < Cd~^ for all multi-index /3, 
then for any if)h G Sh(B 3 d(D)) there exists rjh £ S®(B 3 d{D)) such that 

IIWVi - r lh\\H k {B3 d (D)) < Kh 1 k d 1 \\lph\\L 2 (B 3d (D)), k = 0,1. 
Furthermore, if u = 1 on Bd{D), then r/h = iph on D and 

\\uiph ~ Vh\\H^(B 3d (D)) < A'/i 1_fc d _1 ||^/ l || i 2( B3d(D )\ D ), & = 0,1. 
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For example, rjh = Iih(u>iph) has these properties. 


(P3) Regularized Delta function: 

For any x$ £ f l j, there exists a function S Xo £ G 3 (C2) with support in rj 1 such 
that 


Xh(x o) = / XhS Xo dx, V Xh £ Sh, 


\\5 x J wl , P <Kh 1 iV ( 1 1/,p -> for 1 < p < oo, 1 = 0, 1,2,3. 

(P4) Discrete Delta function 

Let S xo denote the Dirac Delta function centered at Xq, i.e. f n 6 X0 (y)<p(y)dy = 
<p(xo) for any ip £ C(S2). The discrete Delta function PhS X0 satisfies that 


PhS Xo {x) < Kh N e 


I x — x 0 I 


\/x,xq £ SI. 


The properties (P0)-(P4) hold for any quasi-uniform partition with those stan¬ 
dard finite element spaces and also, have been used in many previous works such 
as [23 ESS EH 02. The proof can be found in the appendix of [41] . 

For an element r/ 1 and a point Xq we let G(t, x, xo) be the Green’s function 
of the parabolic equation, defined by 


(2.5) dtG(t, -,xo) + AG(t,-,xo) = 0 for t > 0 with G(Q,x,xo) = 6 Xo (x), 


The regularized Green’s function r(f, a;, xo) is defined by 

(2.6) d t T(-, -,xo) + Ar(-, -,xo) = 0 for t > 0 with T(0,-,x 0 ) = 5 Xo , 


where 6 Xo is given in (P2), and the discrete Green’s function r^-, -,Xq) is defined 

by 

(2.7) d t T h (-,-,x 0 ) + A h T h (-,-,x 0 ) = 0 for t > 0 with r h (0,-, x 0 ) = PhS Xo - 

The functions G(t,x,x o) and Th(t,x,xo) are symmetric with respect to x and Xq- 
By the fundamental estimates of parabolic equations, there exists a positive 
constant C such that m, Theorem 1.6; note that the Green’s function in the 
domain Q is less than the Green’s function in M. N ) 

(2.8) |G(t,x,y)|<C'(t 1 / 2 + |x- 2 /|)- Ar e-^. 

By estimating T(t,x,Xo) = J n G(t,x,y)S Xo (y)dy, it is easy to see that (12.81) also 
holds when G is replaced by T and when max(t 1 / 2 , \x — y |) > 2 h. 


2.3. Decomposition of the domain D x (0,T). Here we present some further 
notations on a dyadic decomposition of the domain Six (0, T), which were introduced 
in [39] and also used in many other articles mmmm- Let R 0 be the smallest 
distance between a corner and a closed face which does not contained this corner. 

For the given polygon/polyhedron SI, there exists a positive constant Kq > 
max(l,f?o) (which depends on the interior angle of the edges/corners of SI) such 
that 

(1) if zq is a point in the interior of SI and B p {zq) intersects a face of SI, then 
Bp(zo) C B 2p (zi) for some z\ which is on a face of SI; 

(2) if z\ is on a face of SI and B p {z\) intersects another face, then B p (z\) C 
BpK 0 { z 2 ) for some z 2 which is on an edge of SI; 
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(3) if Z 2 is on an edge of and H p (z 2 ) intersects another face which does not 
contain this edge, then B p (z2 ) C B p k 0 (z3 ) for some z 3 which is a corner of fl. 

For any integer j > 1, we define dj = 2 ~^~ 3 RqKq 2 . For a given x 0 £ fl, we let J* 
be an integer satisfying dj, = 2~ j *~ 3 RqKq 2 = C*h withC* > max(10, 10k, RqKq 2 /8) 
to be determined later. Thus, J* = log 2 [i?o^o" /(8C*/i)] < log 2 (2+l /h) and J* > 1 
when h < RqKq 2 /(16C*). Let 


Q*(^o) - 

= {(s,t) e 

Qj 1 : 

max(|x — xq|, i 1 ' 2 

) <dj.} 

^*(£ 0 ) — 

= {x S S2 : 

k- 

sol 

VI 


Qj(x 0 ) = 

= {(s,t) e 

: 

dj 

< max(|x — xo 

l,f 1/2 )< 

(x 0 ) — 

= {x £ S2 : 

^ < 

\x 

So | ^ ddj }, 


^j(^o) = 

= {x £ SI : 

|x - 

Sol 

— 2dj} 



for j > 1; see Figure [T] 



Figure 1. Illustration of the subdomains Qj, ttj and Dj. 

For j = 0 we define Qo(xo) = Qt\Qi(so) and flo(so) = £2\^i(so)i and for j < 0 
we simplify define Qj(x 0 ) = flj(x 0 ) = 0. For all j > 1 we define 

Qj(xo) — i(so) U flj(xo) U (xo), 

fl"(x 0 ) = fb,_ 2 (x 0 ) U Slj(xo) U fL, +2 (x 0 ), 

flj'(xo) = flj- 2 (x 0 ) U flj(x 0 ) U fb, +2 (x 0 ), 

Qj(x 0 ) Qj— i(so) Id Qj(x 0 ) U Qj- li(xo), 

Q"(s 0 ) = Qj-2(x 0 ) U Qj(x 0 ) U Q j+ 2 (xo), 

Qj'(x 0 ) = Qj-2(xo) U Qj(x 0 ) U Q j+ 2 (xo), 

Dj(x 0) = U Dj(x 0), 

D"(x 0) = D j _ 2 (x 0 ) U Dj(x 0), 

Dj'(xo) = Dj_s(xo) UT"(xo). 

Then we have 

x, J* 

fir = [J Qj(xo) U Q*(xo) and fl = [J flj(xo) U fl*(xo), 
l=o 1=0 
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We refer to Q*(xo) as the “innermost” set. We shall write j when the innermost 
set is included and ]Tb when it is not. When xq is fixed, if there is no ambiguity, we 
simply write Qj = Qj{x 0 ), Q' = Q'-(xo), Q" = Qj{x 0 ), fij = Clj(x 0 ), fi'- = fi'-(x o) 
and fi" = Q" (x 0 )- 

In the rest of this paper, we denote by C a generic positive constant, which will 
be independent of h, Xo, and the undetermined constant C* until it is determined 
at the end of section m 

2.4. Proof of Theorem 11.11 The keys to the proof of Theorem 11.11 are several 
more precise estimates of the Green’s functions. Let F(t) = r hit) — r(t). Then for 
any Xq £ 12, we have 

{E h {t)v h )(x 0 ) = (F(t),Vh) + (T(t),v h ) 

(2.9) = [ (d t F(s),v h )ds + (F(0),Vh) + {T(t),Vh), 

Jo 

(td t E h (t)v h )( x 0 ) = (td t F(t),v h ) + (td t r(t),v h ) 

(2.10) = [ (sd ss F(s) + d s F(s),v h )ds + (td t T(t),v h ), 

Jo 

with ||F(0)|| i i = ||4 0 — PhbxQ IIl 1 < C (according to (P3) and (P4)). Moreover, by 
the analyticity of the continuous parabolic semigroup on L 1 (f2), we have 

l|r(t)IUi + f|IW)|| L i < q|r(o)|| L i = c\\s X0 \\ L i < c. 

We present some estimates of these Green’s functions in the following lemma. 
The proof of the lemma is the major work of this paper and will be given in the 
next two sections. 

Lemma 2.2. Under the assumptions of Theorem \l.l[ we have 

(2.11) f f (\d t F(t,x,x 0 )\ + \tdttF(t,x,x 0 )\)dxdt < C, 

Jo Jn 

(2.12) \Vd t G(t,x,xo)\ < Gmax(t 1 / 2 , \x — aio|) 3 ^ f or (x, t) £ fl x (0,1). 

The estimates in Lemma 12.21 were proved in [39] for parabolic equations with 
the Neumann boundary condition and in E! for the Dirichlet boundary condition. 
However, their proofs are only valid for smooth coefficients and smooth domains (as 
clearly mentioned in their papers). Later, these estimates were proved in [25] for 
parabolic equations in smooth domains of arbitrary dimensions under the Neumann 
boundary condition with Lipschitz continuous coefficients. Here we are concerned 
with the problem in a convex polyhedron in two or three dimensional spaces under 
the Dirichlet boundary condition with a^- € W 1 ' N+a . 

Proof of Theorem \l.l\ Firstly, from (12.9I) - (12.10I) we see that (11.91) is a consequence 
of (imp . 

Secondly, we can view Eh(t) as an analytic semigroup on L q (UL), defined by 

(E h (t)v)(x 0 ) := / T h (t,x,x 0 )v(x)dx, 

Jn 

whose generator is AhPh- From [49, Theorem 4.2] and [50] Lemma 4.c] (with a 
duality argument for the case q > 2) we know that the maximal L p regularity 
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(fTTTol) 

(2.13) 


holds if the following maximal ergodic estimate holds: 


sup- [ |£7ft.(s)| vds 
t >o t Jo 


<c\\v\\ Lq , VvcL^n), 


LI 


where 

(|S/ l (s)|u)(a;o) := / \T h (t,x,x 0 )\v(x)dx. 

J n 

Let Gtr{t,x,Xo ) be a truncated Green’s function which is symmetric with respect 
to x and xo and satisfies Gt r (t, x, xo) = G(t,x,xo) when ( x,t ) is outside Q*(xo) 
(see U32 Section 4.2] on its construction). Then we have (assuming that Tq is the 
triangle/tetrahedron which contains xo) 


[Ox(0,oo)]\(p»(a:o) 


[Ox (0,1)]\Q* (xq) 


\d t T(t,x,x 0 ) - d t Gtr(t,x,x 0 )\dxdt 


d t G(t , x , y)5 Xo ( y)dy - d t G(t , x, x 0 ) 


dxdt 


|<9 t r(f,x,Xo) — d t G(t, x,Xo)|dxdt 


< Ch 


Qx (l,oo) 


sup \V iy dtG(t,x,y)\dxdt 
[nx(0,l)]\Q,(rc o ) yerft 


/ OO 

t~ 1 (\\T(t/2, ■, ®o)||li + ||G(f/2, -,x 0 )\\ L i)dt [by semigroup estimate] 

= Ch.'S'' ff sup | V y dt.G(t, x, y)\dxdt + C f t~ 1 ~ N / 2 dt [see ()2.8D ] 

■ J JQj(x 0 ) (y,t)&Q'Ax) Ji 


<cJ2^- + c [see mnai 


< c. 


By using energy estimates, it is easy to see 


(|9tr(i,x,xo)| + \d t Gtr(t,x,xo)\)dxdt 


J J Q*{x o) 

< 4 /2+1 (||9 t r(, •,a;o)||L 2 (n x (o, 1 )) + ||<9tG tr (-, •, £o)||l 2 (Qx(o,i))) < C*^ 2+x , 

where the constant G* will be determined at the end of Section 4. Then (12.111) and 
the last two inequalities imply 



\d t r h (t,x,x 0 ) 


d t Gt r {t,x,xo)\dxdt < C. 


In other words, the symmetric kernel A'(x, y) := / 0 °° 9 t T^(f, x, y) — d t G* T (t , x, y)|di 
satisfies 

sup / K(x, y)dx + sup / K(x,y)dy<C, 
yen Jo xeoJn 

and therefore, Schur’s lemma implies that the corresponding operator Mk , defined 
by Mkv(x) = f n K(x, y)v(y)dy, is bounded on L q (VL) for all 1 < q < oo. Let 
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E* t ( t)v(x) = J Q G* T (t,x,y)v(y)dy and note that E* I (t)v(x) < E(t)\v\(x) (because 
G* r (f,x,y) < G{t,x,y)). We have 


sup(\E h (t)P h \v)(x) 
t> o 

< sup(|£ l / 1 (t)P? l |i; - E* r (t)v)(x) + sup (E* r (t)v)(x) 

t> 0 i>0 


< sup (I E h (t)P h - E* r (t)\ H)(x) + sup (E* T (t)v)(x) 

t >0 £>0 


= sup 
£>0 


(IWM) 


f [ I w* 

Jo J n 


,y) 


+ sup(E* r (t)\v\){x) 

t> 0 


d t G* tT {s,x,y) | |-u(y)|di/ds 


< (\PhS x \, |w|) + (M k \v\)(x) + sup(P(t)|i>|)(x) 

£>0 


where 

II supP(t)|u| \\li < C q \\v\\ L i, V 1 < q < oo, 
t>o 

is a simple consequence of the Gaussian estimate (12.81) (Corollary 2.1.12 and Theo¬ 
rem 2.1.6 of BED- This proves a stronger estimate than (12.131) . The proof of (11.101) 
is completed. 

Finally, OD-lfOl imply that the error = PhU — Uh satisfies the equation 
(when u° = = 0) 

(2.14) d t {A~ l e h ) + A h (A^ 1 eh) = PhU - R h u. 

By applying (11.101) to the equation above, we obtain 


(2.15) \\&h || L p ((0,T);L 9 ) < C p , q \\PhU — Rhu\\ L P((0,T)-,Li) < C'p,qh\\u\\ L p(( 0t T)-,W 1 ’i) 

for 1 < p < oo and 2 < q < oo, where we have used the inequality \\PhU — RhuWli < 
G 9 /i||u||t^i, 5 , which only holds for 2 < q < oo in convex polygons/polyhedra. Then, 
by using an inverse inequality and (12.41) . we have 

ll^/iIIL p ((0,T);W 1,9 ) < Ch 1 ||e? l || i p((o ! T);L 9 ) 

< C P}q \\u\\ L p^ 0 T y w i, q ) 

which implies (11.131) for the case 1 < p < oo and 2 < q < oo. 

In the case 1 < p < oo and 1 < q < 2, we define g = \7A~ 1 P/ l f and express the 
solution of <o by (when it° = 0) 

rt 

Vu h = C h g:= / VA~ 1/2 A h E h (t - s)A~ 1/2 V ■ g{s)ds. 

Jo 

In order to prove the boundedness of the operator Ch on L p {{ 0, T); ( L q ) N ), we only 
need to prove the boundedness of its dual operator C h on L p ((0,T); ( L q ) N ). It is 
easy to see that 

{£hg,y)dt = J (§> VA~ 1/2 A h E h (t - s)A~ 1/2 V ■ f}(t)dtjds, 

which gives 

r T 

£ h fj:= / VA~ 1/2 A h E h (t - s)A~ 1/2 V ■ fj(s)ds. 
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If we define the backward finite element problem 

/ 2 lg \ f ~{9tWh,Vh) + (aVw h , Wv h ) = (V • if,v h ), Wv h eS h , 

\ w h {T) = 0, 

then C' h ff = Vi Vh- By a time reversal we obtain, as shown in the last paragraph, 
||Vw/ l || i p'((o jr ); i? ') < Cp,(j||V • vWlp 1 Ill'll Lp' ((0,T);L^ / ) 5 

for 1 < p' < oo and 2 < qf < oo, which implies the boundedness of C' h on 
L p ((0, T); ( L q ) N ). By duality, we derive the boundedness of Ch on L p {{ 0, T); ( L q ) N ) 
and therefore, 

||Vu/i||lp(( 0,T);L9) < Cp,g||fflUp((0,T);L9) 

< C Pt q\\Phf\\LP(0,T)-,W- 1 ’i) 

This proves (11.131) in the case 1 < p < oo and 1 < q <2. 

The proof of Theorem 11.11 is completed (based on Lemma 12.21) . I 

Remark 2.1 In the proof of (11.131) . we have used an L q error estimate of the Ritz 
projection for 2 < q < oo, which can be proved in the same way as used in [36] 
Corollary] by using the kL 1 ’ 9 -stability of the Ritz projection. This fR 1,9 -stability is 
based on an interpolation between these two cases q = 2 and q = oo. The case q = 2 
is trivial and the case q = oo was studied by several authors, such as [36j for r = 1 
and 2D convex polygons (which requires H 2 regularity of the elliptic problem), j37j 
for r > 2 and 2D arbitrary polygons (as a consequence of the L°° stability proved 
therein, which only requires H 3 ' 2+e regularity of the elliptic problem), and [19] 
for r > 1 in 3D convex polyhedra (which requires H 2 and C 1+a regularity of the 
elliptic problem). These essential properties used by [TS] ESI EZ] are all possessed 
by the elliptic problem when the domain is convex polygonal/polyhedral and the 
coefficients are W 1 ' N+a . 

In the rest of this paper, we focus on the proof of Lemma POl 

3. Preliminary analysis 
I n this section, we present two propositions. 

3.1. Superapproximation of smoothly truncated finite element functions. 

In this subsection, we prove the following proposition, which is needed in proving 
Lemma 1^1 

Proposition 3.1. If 0 < lo < 1 is a smooth cut-off function which equals zero in 
fl\D, satisfying \d^to\ < CdT ^ for all multi-index ft such that |/J| = 0, 1, • • • , r + 1 
and d > lOnh, then for any ifh £ Sh there exists Xh £ S^(Bd(D)) such that 

(3.1) d 2 \\R h (uifh) — XhWm + d\\u)i/>h ~ XhWl 2 < Ch\\if> h \\ L 2^ Bd ^y 

Proof. Define 0 < cc < 1 as a smooth cut-off function which is zero outside Bq .8 d{D), 
satisfying that w = 1 on Bo. 7 d(D) and \d l3 u\ < Cd ~^ for |/3| = 0,1, • • • , r + 1. 

First we prove the following inequality 

(3-2) \\un/j h - R h {uj4>h)\\L2 < Chd 1 ||V'/i||i 2 (Bo. 3 d(u > )) 
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by a duality argument. We define </> as the solution of the elliptic PDE 

f —V • (■ aS7(j> ) = v in f2, 

| (f> = 0 on dfl. 

We see that 

(v,u}ip h - R h (uii/j h )) = (dV<t>,'S7(u}ip h - Rh(u}iph))) 

= (aV(0- R h <t>), - Rh(ui/Jh))) 

= (aV(i^-12^),V(w4 -Il^(wV'h))) 

< c||^ - - n h (uijj h )\\ H i 

< Chd~ 1 \\v\\ L 2\\tp h \\ L 2(^B 0 . 3d (D))- 

where we have used the superapproximation property (P2) in section 12.21 and the 
H 1 error estimate: 

H - RhtWm < ChU\\ H * < Ch\\v\\ L 2. 

(13.21) follows these inequalities. 

Secondly, it is noted that the following inequality 

(3.3) \\Rh{u'llJh)\\m(B d (D)\B 0 . 5 d (D)) < C ; d~ 1 \\Rh(ujllJh)\\mB d (D)\B 0 . 3d (D )) 

was proved in Lemma 4.4 of [40] (also see Page 1374 of [39]) as a consequence of 
the discrete elliptic equation 

(aVR h (unj> h ), V?/) = 0, for rj G S° h (B d (D)\D ). 

Let \h = Tlhlu}Rh{toil>h)\ and note that the support of \h is contained in Bo sd(D). 
By using the superapproximation property (P2), we have 

d~ 1 \\Rh(ui/jh) - Xhh 2 + || Rh(uiph) - XhWrn 

< d~ 1 \\R h (uirph) - ujR h (un/j h )\\ L 2 + d~ l \\uR h (wil) h ) - IL h [uR h (uip h )]\\ L a 
+ || R h (uip h ) - uiR h (ujip h )\\ H i + \\uR h (ujiph) ~ Tlh[ujR h {i2jiph)}\\m 

< Cd~ 1 \\Rh{wil>h)\\L 2 (B d (D)\B 0 . 3 d (D)) 

= Cd~ 1 \\R h {uj'tph) - u'i/)h\\L 2 (B d (D)\B 0 . 3 d (D)) (because w = 0on B d (D)\B 0 , 3d (D )) 

(3.4) 

< Chd~ 2 \\ilj h \\ L 2 ( Bo 3d (o)) (as a consequence of (E3) 

and from (13.21) we see that 

II Ulph - Xh ||L 2 < \\ui>h ~ Rh(^h)\\L 2 + || R h (u]l/j h ) - Xh\\L 2 

(3.5) < C'M _1 ||^|| i 2 (Bo3d ( £ j)). 

(13.11) follows immediately and the proof of the Proposition 13.II is completed. I 

Remark 3.1 In the proof of Proposition 13.II we have assumed that d > 10«:d and 
used Bo. 3 d(D), Bo. 7 d(D), Bo.$ d (D) to make sure that their radius differ from 
each other by at least nh so that the superapproximation property (P2) can be 
used. 
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3.2. Local error estimate. The following proposition is concerned with a local 
energy error estimate of parabolic equations. 

Proposition 3.2. Suppose that (j>,dt(j> £ L 2 ((0, T); Hq) and fih £ i? 1 ((0, T); Sh), 
and e = <f>h — </> satisfies the equation 

(3.6) (e t , X ) + (aVe,V X ) = 0, V y £ S h , t > 0, 

with 0(0) = 0 in Qj" ■ Then for any m > 0, there exists a constant C m > 0, 
independent of h and dj, such that 

l|et||i 2 (Q J ) +dj 1 |l^ 7e IU 2 (Q J ) 

(3.7) < 0)) + *i(n Kt -</>) + Hj (e) + df\\e\\ L * W p), 

where 

0)) = ^7 1 ||<A/ l (0)||L2(ny') + II <t>h (0)11 jthjj'")) 

(IIh0-0) = dj||V9t(IIfe0 — 0 )||i,2(qw) + ||9 t (II^0 - 0 )||l2 ( q'") 

+ d7 1 ||V(n /l 0 - 0 )||l2(q///) + d” 2 ||n^0 - 0 ||l2(q"'), 

Hfie) = {h/dj) m (\\e t \\ L2wn +dJ 1 \\Ve\\ L ^ QT) ). 

Before we prove Proposition 13.21 we present a local energy estimate for finite 
element solutions of parabolic equations. 

Lemma 3.3. Suppose that 4>h{t) £ Sh satisfies 

ipt^hi X) + (aV0fc, Vx) = 0, for X £ S%(Q"), t £ (0, d 2 ], 

(9*0/!, x) + (aV0/ l , Vx) = 0, /or x £ t £ (d 2 /4,4d 2 ). 

XTien for any m > 0 f/iere exists C m > 0, independent of h and dj, such that 
||5t0/,||z,2(Q ;j ) + d7 1 ||V0 fc || i 2 ( g j) 

< C m (||V0/,(O)||i2(n") + d“ 1 |l0ii(O)IU 2 (n")) 

( i \ m 

(119*0/! ||l 2(Q») + d7 1 ||V0/,|| L 2(Q//)) +C' m d7 2 ||0/ l || L 2(Q//). 

Proof. Note that Qj = [1L,- x (0,d 2 )] U [Dj x (d 2 ,4d 2 )]. We first present estimates 
in the domain Qj x (0, d 2 ) and then present estimates in the domain Dj x (d 2 ,4d 2 ). 

Let a; be a smooth cut-off function which equals 1 in Qj and vanishes outside Qj , 
and let ui be a smooth cut-off function which equals 1 in Q'" and vanishes outside 
Qj" with 

(3.9) \d^uj\<Cd~ m and \d^\<Cd~ W . 

Let v h = II hiuxfh) £ S%(Q"'") so that v h = <f> h in Q'J, satisfying (due to the 
superapproximation property (P2)) 

IKIIz, 2 < C'||0/t|U 2 (fi""), 

||Vu/ l || i 2 < C , ||V0 fc || i 3 (nr) + CdJ 1 ||0/!||L2(Q////) 

and 

(d t v h , x) + (aVv h , Vx) =0, V X £ S°(ft") • 
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It follows that 

i-^lla^H 2 + (u} 2 aS7v h ,S7v h ) 

= [( d t v h ,uj 2 v h ~Xh) + (a\7v h ,V(R h (u] 2 v h ) - Xh))] 

+ [( uj 2 aX 7 vh,Vvh ) - (aVvh,V(u} 2 v h ))] 

< [CWdtVhW^WvhW^hdJ 1 + C\\\7v h \\ L 2\\v h \\ L 2hd~ 2 ] + C{uaMv h , 2v h Vu ), 
where we have used (P2) and Proposition 13.11 and from (13.91) we see that 
(uaVwft, 2u^Vw) < (\u>aVv h \, 2 1-y*,|)<i~ 1 < C\\ioaS7v h \\ L 2\\v h \\ L 2d~ 1 . 

The last two inequalities imply 

ll < / , fi||i“ I ((0.d|)xL 2 (n i )) + ll^ 7 ^llL 2 ((0,d 2 );L 2 (nj)) 

< (711^(0)11^2(0////) + C'||9 t ^|| i 2(( 0id 2). i 2(O'"'))/l 

(3.10) + C\\ V0? l || i 2(( Oid 2). i 2(o'"'))/ld J - 1 + C'||^f i || i 2(( 0jd 2) ;L 2(o'"'))<i :) - • 

By using Proposition 13.11 again, we derive that 

II u 2 d t v h \\ 2 L 2 + ^^(u 4 aVv h ,Vv h ) 

= [( d t v h ,u] 4 d t v h - Xh) + (aVv h ,V[Rh(oj 4 d t v h ) - Xh])] + (4w 3 aV^, d t v h Vu) 

< CWdtVhW^hdJ 1 + C\\Vv h \\ L 2\\dtVh\\L*hd- 2 + CWuVvhWvWuPdtVhWvdr 1 

< C\\dtVh\\L2hdj 1 + C\\S7vh\\L2dj 2 + ^||w 2 9tu/ 1 ||| 2 , 
which reduces to 

||w 2 9 t U/ l || 2 2(( 0id 2); i 2(O)) 

< C'||Vu^(0)|||2(o) + C'||3 t ?;/ ( || 2 2((o j d2);i, 2 (o))^j 1 + ^ll^/Jl 2 ((o,d 2 );L 2 (n))^j 2 - 

The inequality above further implies 

l|5t</>/ l ||L2(( 0d 2) ;i 2(O J .)) < C(||V</>/i(0)|| i 2(O/'") +d j 1 ||<(>?l(0)|| L 2(o'/'/)) 

+ C\\d t <t>h\\L 2 ((o,d 2 j );L 2 (n' j , "))h 1 ^ 2 dj ^ 

(3-11) + Cd j 1 (ll^fi IIL°°((0,rf 2 )xL 2 (Q' /,/ )) + l|V</>/, || L 2 ((o,,j 2 );.L 2 (n"")))- 

With an obvious change of indices (from f V" to Q!" on the right-hand side, and 
from Q to ST' on the left-hand side), (13.10l) - (13.11l) imply 

ll^llL“(( 0 ,d 2 );L 2 (n')) + IIW/i||.L 2 (( 0 ,d 2 );Z, 2 (n')) 

< (711^(0)11^2(0'") + Ch\\d t (t>h\\L 2 ((o,d 2 j )-,L 2 (n l j ")) 

(3.12) +Chdj 1 ||V</>/i|| i 2((o d 2) ;i 2(o///)) +Cdj 1 II z, 2 ((o,d|) ; z,2(o'")) - 

and 

ll^t<( ) ?i||L 2 ((0,d 2 );L 2 (n j )) < C(ll Wh(0)llL 2 (n'.) +dj 1 ||<^/i(0) 11^2(0')) 

+ Ch 1 / 2 d J ! ||5t^||i2(( 0(i 2). i 2(O')) 

+ (ll‘MlL“((0,dJ);.L 2 (fi')) + II W/i||,L 2 ((0,d 2 );L 2 (n')))- 


(3.13) 
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In the same way as we derive (13.12|l - (13.13ll . by choosing w(x,t) = ui(x)uj 2 (t) 
with uji = 1 in £>'■, = 0 outside D", W 2 = 1 for t G (d|,4d|) and Il >2 = 0 for 

t G (0,d^/2), we can derive that 

l|w^7li“((C),4d 2 );L 2 (n)) + ||wV^/,|| i 2 (( 0 i 4 d 2 J ;i 2 (f 1 )) 

< C'/l|| 9 t 0 h|| i 2 (( d 2 / 4 )4d 2 ) ;i 2 ( £ )///)) 

+ Chdj 1 ||V^|| i 2 (( d 2 / 4 j 4 (i 2 ) ;L 2 ( £ )///)) 

(3-14) + 1 ||^Ii||L 2 ((d2/4 > 4d2);L 2 (D'")) 

and 

ll | 5t^ft||L 2 ((d 2 ,4d 2 );L 2 (D 3 )) < Chf^dj ^ ||9t^fc||z,2((d2/ 4l 4d2). i 2(£,/»)) 

(3.15) + C'dT 1 (||G%|| i 

“((0,4d 2 );L 2 (n)) + ll w W?i|l,L 2 ((0,4d 2 );.L 2 (fi))). 

By noting the definition of w and u;, we have 

\\dt<f>h\\L 2 (Qj) < C'(l|5t < /'Ii||L 2 ((0,d 2 );L 2 (n 3 -)) + || L 2 ((d 2 ,4d 2 );L 2 (_Dj))) 

l|V</>/i||i 2 (Qj) < C(\\u}\/(/)h\\L 2 ((0,d 2 j )-L 2 (n)) + ll a ''^ 7< /’/i||L 2 ((0,4d 2 );L 2 (n))) • 

With the last two inequalities, combining (13.10[l - (13.15l) gives 

\\dt4’h\\L 2 {Q j ) + dj 1 ||V(/)/ l || L 2(Q i ) < tf(||V^(0)|| L 2 (n „,) +dj 1 ||^(0)|| L 2 (n ///)) 

+ c {\\dt<t>h\\L*{Q'p + dj 1 \\S/4>h\\L 2 (Q'J')) 

+ C^7 2 || ( / > h||L 2 (Q'"). 

Iterating the inequality above and changing the indices, we derive (13.81) . I 

Now we are ready to prove Proposition 13.21 

Proof of Proposition [X2I Let u)(x, t) be a smooth cut-off function which equals 1 
in Q" and vanishes outside Q ”', and let f> = ui<j>. Then (f>(0) = 0 and we have 

{d t {4>-cfh),x) + (aV(0- <f h ), Vx) = 0, for x e S^ft'), t € (0,d%], 

(/>h),x) + (aV(0- (j>h), Vx) = 0, for x e S ° h (£>'■), t G (d|/4,4d|). 

Let (fh G 5h be the solution of 

(3.16) [d t {4>- 4>h),Xh) + (aV(0- 0h), Vxh) = 0, for xh 6 S’/, 
with ^(0) = 11^(0) = 0 so that 

(3 - 17 L 

(5 t (</>/! - <fh),Xh) + (aV(^ h - </>/>), Vxh) = 0, V x/, 6 ^(Sl'), te(0,d-], 

(3-18)_ 

- <t> h ),Xh) + (aV(^ h - 4> h ), ^Xh) = 0, V x/i e 5^(£>'•), t G (d|/4,4d|). 
Substituting x/i = Ph4> — 4>h into (13.161) we obtain 

j|^-^||7=o ((0iT);i 2) + [ (aV{<j>-(i> h ),V((j>-(j) h ))dt 


0 
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[ {d t {4>-4>h),4>-Ph4>)dt + f (aV(0 - 0/0, V(0 - P h (/>)) dt, 
Jo Jo 


which implies 


110 _ 0/Jl°“((O,T);L 2 (Q t ) + IY(0 — 0h)ll L 2 (Q t ) 

(3.19) < C\\dS-M\\ LHQT) U\\ L2{QT) + C||V^||| 2(Qt) . 

Substituting \h = d t (Ph4> — 0/0 into (13.161) we obtain 


||<9 t (0 — 0 /i)||l2(q t ) 

r T 


sup (aV(0 — 0/0, V(0 — 0/0) 

0 <t<T 


= (<9 t (0-0/0, <9/(0--P/i0))df + / (aV(0 — 0/0, Vd t ((j) — P/i0))df 

Jo Jo 

< ||<9/(0 - 4>h)\\L 2 (QT)\\dt4>\\L 2 {Q T ) + II V(0 — 0/ 1 )||L 2 (Q T )||Va t ^|| L 2 ( Q T ), 


which implies 
(3.20) 

\\dS-M\\l HQT) < G||St0||| 2(QT) + C||V(0 - 4 h)\\l*(q t )\\VM\\l>(q t ), 
It follows that 


ll^t(0 — 0h)lli 2 (Q T ) + dj 2 ||V(0 — 0/i)||i2(Q T ) + d■ 2 \\(j) - 0fc||i«.(( O)T ) ;i 2) 

^ 2ll^t(0 — 0^)lli 2 (Q T ) + Cdj 1 ||0Hl 2 (q t ) + Cdj "|| V0|||2(q t ), 

+ C'||9t0|||2(Q T ) + -dj 2 1|V(0 — 0/OIIl 2 (q t ) + C^||Va t 0|||2 (Qr) . 
which in turn produces 

||<9 t (0 - 0/011 L 2 (Qt) + dj 1 ||V(</> - 0/OH L 2 (Qt) + dj 1 10 - 0/i ||l oo ((o,t) ; l 2 ) 

< Cd~ 2 U\\ L . {QT) + Cd- 1 ||V0|| i 2 (c3T) + C||a t ?|| i 2 (QT) + Cd J \\Vd t ^\\ L 2 {QT) . 

By applying Lemma [3731 to (13.17l) - (|3.18l) and using the inequality above, we derive 
that 


II dt{<t>h - 0/0II z. 2 (Q j ) + dj 1 ||V(0? l - 0/ 1 )||l2(q 3 .) 

< C m ( |V(0h - 0h)(O)||x,2( n //) + d^W^h - 0h)(O)|| L 2 (n //)j 

( i \ 771 _ 

(II 9t{<j>h - 0/i)||i 2 (Q") + dj 1 IIv(0ft, - 0/i)H L 2 (Q'/)) 

+ C m dj "H0/! — 0/i||l 2 (Q'') 

< C m (||V0fe(O)||i2( n ») + rf7 1 ||0/i(O)||L 2 (n'')) 

( h \ m 

(ll^e|U 2 (QC) + rf 7 1 |Ye|| L 2 ( Q// ) ) + C' m d“ 2 ||e|| L 2 ( Q// ) 

( i \ 171 _ _ _ _ 

(11^(0 — 0/O I|l 2 (<2") + dj 1 1| V(0 - 0/O||l 2 (Q")) 

+ C m e^ 1 1| 0 — 0 /i||l“L 2 (Q") 

< (||V0/ l (O)|| i 2 (ny) +d- 1 ||0/ l (O)|| i 2 (n ,, ) ) 
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( h \ m 

-j-J (ll^ e lk 2 (Q'/) + dj 1 ||Ve|| £ 2 (Q//)) +C m dj 2 ||e|| £ 2(Q») 

+ Cmdj 2 \\Hl 2 (Q t ) + Cmdj 1 ||V0|| L 2(q t) + Cm\\dt4>\\L 2 {Q T ) 

+ C m dj\\Vdt4\\vnQ T ) 

The last two inequalities imply 
(3.21) 

II e tIIz, 2 (Qj) + dj 1 1|Ve||i 2 (g j ) < C m (I,(MO)) + Xj(cf)) + Hj(e) + dj 2 l|e||z, 2 (Q'")) ■ 

We have proved that any e = fa — (/> satisfying (13.61) also satisfies (13.211) . Since 
e = fa — II h,(j> — {<j> — II h4>) and </>(0) — II/j^O) = 0 in Q'”, we can replace fa by 
fa — Hh<fi and (j) by <j> — Hh(f in (13.211) . Then (13.71) follows immediately. I 

4. Proof of Lemma 12721 
Now we turn back to the proof of Lemma 12.21 

4.1. The proof of (12.121) . In this subsection, we present several local energy 
estimates for the Green’s function, the regularized Green’s function and the discrete 
Green’s function, which then are used to prove (12.121) . These energy estimates will 
also be used to prove (12.111) in the next subsection. In this subsection we let T = 1 
and fix xo £ f l. We write G and T as abbreviations for the functions G(-,-,x o) 
and T{-,-,x 0 ), respectively, when there is no ambiguity. We use the decomposition 
of Section [?~3l for all j > 1 (not restricted to j < J*) and so we do not require 
h < R 0 Kq 2 /(16C*) in this subsection. 

Lemma 4.1. For the Green’s functions T, G and T ? t defined in m - \2. 7| ), we have 
the following estimates: 

(4.1) 

2 

£ df +k - 1+N/2 (\\V k d l t G(-,-,x 0 )\\ LHQj{xo)) + ||V fc 5|T(-,-,a;o)|| I ,2 ( Q.( a . o)) ) < C, 

l,k —0 

(4.2) 

||V 2 G(-,-,aio)|| L <x ) ,2(u fc <.Q fc ( a;o )) < CdJ N/2 ~\ 

(4.3) 

||V Xo atG(v,*o)||£- W , ( x 0 )) < cd- N ~ 3 , 

II d t G(-, ■, ico)||Li(nx(i,oo)) + || td tt G(-, uZcOlUqnxfi.oo)) 

(4.4) 

+ \\d t T{-, -, £o)||z,i(nx(i,oo)) + l|^ tt T(-, •, xo)||Li(nx(i,oo)) < C, 

Proof. For the given Xq and j, we define a coordinate transformation x — Xq = djX 
and t = d?t, and define G{t,x) := G(t,x,x o), a(x) := a(x). Via the coordinates 
transformation, we assume that the sets Qj , Ql, Qj, Ql and f 1 are transformed 

to Qj, Q'j, Qj, Qj and Q, respectively. Let 0 < un{x,t) < 1, i = 0,1,2,3, be 
smooth cut-off functions which vanishes outside Ql and equals 1 in Qj. Moreover, 
uii equals 1 at the points where a3j + i ^ 0, and |Vu5j| < C, \d t 0 Ji\ < C for i = 0,1,2, 3. 
Since U k>jQ'k U Q* is of unit size, there exists a convex domain D = B p (z) D Q D 
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Ufc>, 0' fc U fl*, with 4 < p < 8 Kq, which belongs to one of the following cases (there 
are only a finite number of shapes for D ): 

(i) z £ Cl, p = 4, and B^(z) has no intersection with the boundary of Cl, thus 
B p (z) nfi = B p (z), 

(ii) z is on a face of SI, p = 8 and B$(z) has no intersection with other faces of 
Cl, thus B p (z ) fl Cl is a half ball, 

(iii) z is on an edge of Cl, p = 8 Kq and B$k 0 ( z ) h as n0 intersection with any 
closed faces of Cl which do not contain this edge, thus B p (z) D O is the intersection 
of a ball with a sector spanned by the edge, 

(iv) z is a corner of Cl and p = 8Kq < Rq, and B p (z) fl Cl coincides with the 
intersection of the ball B p (z) with the cone spanned by the corner 2 . 

Note that D x (0,16) contains Q', and consider u>iG, i = 1,2, which are solutions 
of 

(4 - 5) _ 

<%(wiG) — Vs • (aVs(oJiG)) = u} 0 Gdpzi + woGVs • (aVscci) — Vs • (2awoGVsWi) 

and 

(4.6) _ 

dt(u 2 G) - Vs ■ (aVs(w 2 G)) = wiG<%w 2 + wiGVs • (aV 2 w 2 ) - Vs • (2dwiGVsw 2 ) 

in the domain D x (0,16), respectively, both with zero boundary/initial conditions. 
Since D is a convex domain, for p = N+a and po = Np/(N+p) so that W 1 ,P (D) c —> 
L°°(D) and W 1 ,P 0 (D) c —>• L P (D), the standard L p ((0,16); W 1 ,P (D)) estimate of 
(14.51) (the inequality (|2.3[) - (12.4[) with p = q, Lemma 12711) gives 

II Vx( w lG)|| I/P ^o, 1 6);I,i>(5)) 

— GllwoGll^p^Q^gj.^p,,^^ + G||uoGVsa|| iP (( 0il6 ) ;i p 0 (j5)) + G||woG|| i i.(( 0 jl6 );iJ>(5)) 

< G||woG|| l „^ 0i16 ) ; s,!>( 5))(G + G||Vsa|| L N(5)) 

< G||woG|| lp ^ 0 

and the maximal L p regularity of (14.61) yields that (see inequality (12.31) . Lemma 

m 

ll%( w 2G))|| i p(( 0jl6 ) ;i p(5)) + ||Vs • ( a Vs(w2G))|| iP ^ 0jl6 j. iP ^^ 

< G||wiG|| iP ^ 0)16 -j. iP ^^ + GHwiGVsall^p^Q^gj.^p^jj 

+ G||Vs(wiG)|| iP ^ cl;16 ). i p(£ l ^ 

< C llwiGll^p^Q^g^.^po^^^G + G||Vsa|| iP (5)) + G||Vs(wiG)|| i p^ 016 ^. i p^^ 

— H^ 72 ( Wl ^')llLP((0,16);L! , (i5))‘ 

By using m-m, we have 

2 

II Vs(w 2 G)|| i oc ( 5 ) + E l|V|(S 2 G)|| i2(5) < G||Vs • (aV 5 (3 2 G))|| iP(5) . 

k—Q 
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The last three inequalities imply that 

2 

l|Vs( w 2G)|| i 2(( 0>16 ) ;i ao(5p + ^ II^x( W2 ^)IIl 2 ((0,16);L 2 (5)) 

k -0 

+ R(w 2 g)|| L n +“((0,16);L n +“(Z3)) 

< C || w 0 G || iN+ „^ 0 16 ). iJV+a (5)). 

Similarly, replacing G by d^G and d~G in the above estimates, respectively, one 
can derive that 

2 2 

^||V 2 ( W3 5 ? G)|| L2 (( 0il6 ) ;inO ( 5 )) + ^2 ll^( w 3<9jG )|| i2 (( 0 16 ) ;i 2(5)) 

1=0 l , k =0 

< G||woG|| iJV +<»(-( 0 16 ) ;i 2v+o,(5p < C ll w oG|| io0 (( 0j i6); i oo(5))- 
Since W 3 G = 0 at t = 0, it follows that 

II V^c^-(co 3G') II z J °°(( 0 , 16 ) ; i°°(i5)) — G|| ^xdj(uj3G)\\ L 2^ 0ie y Lao ^^ 

— < ^II W ° ( ^IIl“((0,16);L oo (5)))’ 

and 

2 2 

El|V£(S 3 G)|| L 

“((0,16);L 2 (U>)) ^ C E ll^t'^a:( u; 3G)|| i 2( ( -Q 46);L 2 (£))) 
fe=0 fc=0 

— ll W °^llL“((0,16);L“(i3)) ‘ 

Moreover, from the last three inequalities, we have 
2 2 

IIv^IUoo ( q j) + EHv|g|U 2( q 5 . ) + E I$v|g|| £ 2(§j) <c\\g\\ l ^, y 

fc=0 l , k =0 


Transforming back to the (x, t) coordinates, we see from the last two inequalities 
that 


d^dtGW^Q^ + j^d) N/2 \\S7 k G\\ L ~,2 (Qj) + J2 d f +k 1 N/2 \\^ k G\\ L . {Qj ) 

k —0 l,k=0 

< C\\G\\ L ~ Wj) < Cd ~ N , 

where we have used (12.81) in the last inequality. By the symmetry of G with respect 
to x and xo we also get 

d J 3 ||V a:o a t G(-,-,xo)|U=o ( Q. ) < C\\G\\ l <*>(q>.) < CdJ N . 

This proves (14.1I) - (14.8I) for G. 

By using the expression 


(4.7) 


r(x, t ; Xq) = / G(x,t;y)S(y,x 0 )dy. 


one can derive the same estimates for T: 


l|V fc a'r(-,-;a;o)|| z , 2 (Q.) < 


/ HV k d‘G(;-;y)ll L2(Q .)\S(y,xo)\dy 
J n 
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< [ Cd~ 2l ~ k+1N/2 \S(y,xo)\dy < Cd~ 2l ~ k+1 ~ N/2 . 

J n 

Finally, we note that 

l|9 t r(t,-, xo)||l~ +i||d tt r(t, •,x 0 )||l~ 

< C\\d t T(t, -,iEo)||ir 2 + Ct\\d tt T(t, -,x 0 )\\ H 2 

< C\\d t AT(t, ; X 0 )\\ L 2 + Ct\\d tt AT(t, ; X 0 )\\ L 2 [H 2 estimate, Lemma 12.1| 

= C\\dttT(t,-,x 0 )\\ L 2 +Ct\\dtt t r(t,-,x 0 )\\ L 2 

< Cr 2 \\T(t/2,-,x 0 )\\ L 2 [semigroup estimate] 

(4.8) < Ct~ 2 ~ N/ 2 , 

which implies the first part of (14.411 and the second part of (14.411 (the estimates of 
G ) can be proved in the same way. 

The proof of Lemma mi is completed. I 

When max(t 1 / 2 , \x — cco|) < d\ the inequality (12.1211 follows from (14.311 . When 
max(t 1 / 2 , \x — xo|) > di, the estimate ||Va; 0 9 t G(-, •, x 0 )\\L^(Q :j {xo)) ^ C can be 
proved directly (without using the scale transformation) in the same way as above. 

4.2. Proof of (12.1111 . The proof is also based on Lemma 14.11 

First we consider the case h < ho := (RqKq 2 /(16C*) and let T = 1. The basic 
energy estimates of the equations (ESD-(EZD yield 

\\d t Th\\mQ T ) + m\ L 2 (QT) < ||vr„(o)|| L2(n) + ||vr(o)|| i 2 (n) 

= \\VP h 6 X0 \\ L 2( n) + ||V4 0 || L 2 (n) < Ch- X - N l 2 

and 

II WU 2(Qr) + \\d tt T\\ LH Q T) < ||va t r /1 (o)|| L 2 (n) + || w(o)|| L 2 (n) 

= \\VA h P h 6 X0 \\ L 2 (n) + ||VAf X0 || La(n) < Ch ~ 3 ~ N / 2 , 

which imply 

||^f|U 2 (q ,) + \\tdttF\\ L 2 iQt) < ch - 1 -^ 2 + cd 2 h - 3 ~ N ' 2 

< ch- 1 -^ 2 + cc 2 ^ 1 -^ 2 . 

Hence we have 

\\dt-f IIl^Qt) + II^«F|Ui (Q t ) 

< Gd^iWdtFW^Q^ + \\td tt F\\ L 2 (Qt) ) 

+ ^Cd J 1+JV/2 (||a t F|| L 2 (Qj ) + \\tduF\\ L 2 {Qj) ) 

3 

< CCl +N/2 + £ Cd) +N/2 (\\d t F\\ L 2 {Qj) + 11^11^2^.)) 

3 

(4.9) < CCl +N/2 + C7C, 
where 

(4.10) 


K == IYUMoH + \\dtF\\ L 2 (Qj) +d 2 \\d tt F\\ L 2 {Qj) ). 
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We proceed to estimate K,. We set e = F {<f>h = T^ and cj> = T) and e = dtF 
(cj>h = dtTh and cf> = dtT) in ([3.611 (Proposition 13.21) . respectively, and note that 
T(O) = d t r(0) = 0 on f2'-. We obtain that 

d J 1 \\^ /F \\L 2 (Q j ) + \\dtF\\ L 2(Q.) + (PjWdttFW^^Q.) 

(4.11) < C(Ij + Xj + Hj + d■ 2 ||P||l 2 (Q')) 

where 

Ij = ||VP/X 0 ||L=(fi') + dJ 1 \\Ph5 Xo \\ L 2(ci' j ) 

+ d 2 \\V A h P h 5 X0 \\ L 2(n'.) + dj\\A h P h 5 X0 ||i, 2 (n'), 

Xj = dj\\vd t (u h T - r)|| iawj) + \\dt(n h r - r)|| i2(Q ,) + dr^vc^r - r)|| L2WJ) 

+ d- 2 ||n,r - r|| L 2 ( Q/) + - r)|| L2WJ) + d 2 ||5 tt (n^r - r)|| L2(QJ) 

h 3 = (V^nil^FlU^) + d7 1 ||VF|| i2(Q; .) + d 2 j \\d tt F\\ L2{Ql . ) + d j ||W|U 2W ' ) ) . 
By noting the exponential decay of PhS Xo {y) (see (P2) in section [2~2li we derive 
I] < Chd ~ 2 ~ N/2 , 

Xj < (d J /i + /i 2 )||v 2 5 t r|| L2(Q - ) + (dfh + dfh 2 ) ||v 2 r|U 2( Q }0 

+ (d 3 J h + d 2 J h 2 )\\V 2 d tt T\\ LHQ7) 

<ChdJ 2 ~ N ^ 2 , [bv using Lemma |4.1| 

Hj < {h/dj) m (\\d t F\\ L ^ QT) Pd^llVF\\ l2(Qt) +d 2 j \\d tt F\\ LHQT) + dj\\Vd t F\\ L2{QT) ) 

< C'(V^) ro (||^Jx„||iL 1 +dJ 1 \\P h S X0 \\ L , +d 2 j \\A h P h 6 X0 \\ H i +d j \\A h P h 6 X0 \\„) 

< C(/i/d i ) m (/i- 1 - Ar/2 + dj 1 h~ N / 2 + d 2 h~ 3 ~ N ^ 2 + djh~ 2 ~ N /2 ) [by (P3)-(P4)] 
<ChdJ 2 ~ N ^ 2 , [by choosing m = 4 + N/2] 

Therefore, by (14.101) . 

1C < E Chd f + E CdJ 1+N/2 \\F\\ i2(g j) 

J 3 

(4-12) <C + C'^d7 1+JV / 2 ||F|| i2(Q , ) . 


To estimate ||P||l 2 (q'), we apply a duality argument. Let w be the solution of 
the backward parabolic equation 

{ — dtw + Aw = v in f2, 
w = 0 on dfi, 

w(T) = 0 in H, 


where v is a function which is supported on Q' and ||t||^ 2 (q') = 1. Multiplying the 
above equation by F 1 with integration by parts we get 


(4.13) 



Fvdxdt = (P( 0),u>(0)) + 



atj dj Fdiwdxdt , 
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where 

(-F(O), u>(0)) = (PhSx 0 —Sx o ,w(0)) 

= ( PhSx 0 - S Xo ,w( 0) - Ihw( 0)) 

= (PhS Xo ,w(0) - I h w( 0))n» + {PhSx o - 4 OI w(0) - 4w(0))(oy)c 
:= Ii +22- 


By using the exponential decay of P^i5 Xo (see (P4) of section 2) and the local 
approximation property (see (P2) of section 2), we derive that 

N < Ch\\Ph5 Xo ||L 2 ( n'-' ) (dJ 1 11^(0)||z,2 (f!) + ||V«;(0)|| ia(n) ) 

< Ch~ N / 2+1 e~ Cd P h (d j 1 \\v\\ L 2(,n+2)/(,n+4 )(q/) + IM|l 2 (Q')) 

< C7fc-^/ a + 1 e- Cf ^/ fc ||«|U. (a j ) 

< C(dj/h) 1+N ^ 2 e~ Cd P h h 2 dJ 1 ~ N/2 

(4.14) < Ch 2 d ~ l ~ N/2 , 

|2 2 | < C'||4 o || L 2||w( 0) - 4w(0)||i2( (n //)c) 

2 

(4.15) < Ch 2 - N ' 2 ^“ 2 II V fc w(0) || i 2 ((Q po ) . 

fc =0 


To estimate ||V fe u;(0)|| i 2((Q//jo), we let Wy be a set containing (f2") c but its 
distance to 11' is larger than dj /8. Since 

^x w (. x i 0) = f f S7^.G{s,x,y)v(y,s)dyds, 

Jo Jn 

by noting the fact 

\x — y\ + s 1 / 2 > dj/ 8 for x G Wy and (y, s) G Q' 
and using (IT21) . we further derive 


(4.16) 


£d*- 2 ||V fc ttf(->0)IU’ W) 


k—0 


^ C ^l dk j 2sU P\\^ kG (-r,y)h^(\J k<j _ 3 Q k (y))\\v\\ L HQ' 
yen 


k—0 




fc-2 i-JV-fe+JV/2i 




fc =0 

2 




fe-2 i—N—k+N/2iN/2+l 


d^^\\v\\ LHQ/j) =Cd 


-1 

’J ' 


fe =0 


From (I4.14I) - (I4.16I) . we see that the first term on the right-hand side of (14.131) is 
bounded by 

|(P(0),«;(0))| < Ch 2 d/ N/,1 ~ x + Ch 2 - N / 2 d/ 1 < Ch 2 ~ N/2 d~ 1 , 


(4.17) 
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and the rest terms are bounded by 

N 

d t Fwdxdt + y 

i,j =1 


Qi 


Qt 


aijdjFdiwdxdt 


N 


i,j =1 


dtF(w — Hhw)dxdt + // a,ijdjFdi{w — Ilhw)dxdt 


(4.18) 


< yC\\^M\^ m {h 2 \\d t F\\ L , (Qi) + h\\VF\\ L * m ). 


Moreover, to estimate || V 2 w\\l 2 (Q') we consider the expression 

V 2 x w(x,t)= [ I V x G(s — t, x, y)v(y, s)l s>t dyds. 

Jo Jo. 

For i < j — 3 (so that d, > dj ), we see that w(x,t) = 0 for t > 16d 2 (because v is 
supported in Q/); di/2 < \x — y\ < Mi and s — t < df for t < 16d 2 , (x, t) £ Qi and 
{y,s) £ Qj. Therefore, ( x,t ) £ Q'i(y) and we obtain 

ll V2 HI L 2 (Q') < SU P l|V 2 G(-, ■,y)\\L 2 (Q , i (y))\\ v \\L 1 (Q' j ) 

< 2+1 \\v\\ L , m 

< C{d j /d i ) N/2+1 < < =^-. 

For i > j + 3 (so that di < dj), max(|s — t\ 1 / 2 , \x — y\) > dj+i for (x,t) £ Qi and 
therefore, 

II V 2 w|| Z/ 2 (q/) < sup ||V 2 G(-, ■, 2/) 1 u fc <,- +1 Qk(v) IU 2 (Q') IMUdQ') 

yen 

< Cdi sup ||V 2 G(-, -, 2/)l|L°°' 2 (Ufc < 3 . +1 Qfe ( y )) ll 1 ’IU 2 ( Q( ) ^ /2+1 

< Cdid7 N - 2+N/2 d" /2+1 = 

3 3 ^ 

Finally for \i — j\ < 2, applying the standard energy estimate leads to 

IMU 2 ((0,T);ff 2 ) < G||u|| L 2(q t ) = G. 

Combining the three cases, we have 

(4.19) ||V 2 w|| j c 2 (q/) < Cmm(di/dj,dj/di) := Crriij. 

Substituting (|4.17I) - (14.19D into (14.1311 gives the estimate 

(4.20) ||F|| iawJ) < Ch 2 d~ Nl2 ~ 3 + Cym ij {h 2 \\d t F\\ L2(Qi) + h\\ VF|| L2(0i) ), 

which together with (14.121) implies 

/ h \ 2-JV/2 

£ < G + ( - j + cydf 2 - 1 ymij{h 2 \\d t F\\ L2m + /i||VF|| i2(Q<) ) 

j ' 3 ' j *,i 

< G + Cy (h 2 ||<9 t F , || L 2 (Qi) + h\\ VF|| L 2 (Qi) ) y df^mij 
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<C+CJ2 (h 2 \\d t F\\ LHQi) + h\\\7F\\ LHQi) ) d f 72 ” 1 
<C+C (h 2 \\d t F\\ L 2 { Q^ + /i||VF|| i2(Q .)) (CU)^ 2 " 1 
+ Cj24 +N/2 (\\dtF\\ L ^ Qi) +d- 1 \\VF\\ LHQi) ) (j) 

i ' ' 

< c + CC: 1+N/ 2 + cJ24 +N/2 (\\dtF\\ LHQi) + d- 1 ||VF|| i 2 ( g. ) ) fb.) 

i \ 1 / 

<C 2 + C 2 C7 1+N/2 + C 2 C~ l K 

for some positive constant C 2 . By choosing 

(4.21) C* = max(10, 10 k, R 0 Kq 2 /8) + 2C 2 , 

the above inequality shows that K, <C. 

Returning to (14.91) . the boundedness of /C implies 

(4-22) \\dtF\\ L i(Q T) + \\td tt F\\ LHQT) < C. 

From (14.101) we also see that, the boundedness of /C implies 

\\d t F\\ LHQj ) + \\tdttF\\ LHQj) < CdJ 1 ~ N/2 . 

Since fl x (1/4,1) C Udj>i/ 2 Qj, it follows that 

ll^^ 1 |lL 2 (nx(l/4,l)) + l|t^ttF 1 ||L2(Q x ( 1 /4 il )) 

- (II^UI L 2 (Q,-) + II^M^II L 2 (Qj)) 

dj> 1/2 

< Y, Cd~ 2 ~ N < c. 

dj> 1/2 

The above inequality and (14.81) imply 

l|9tr/ l ||i 2 (nx(i/4,i)) + ||tdt t r h || i 2 (n>< ( 1 /4 il )) < c. 

Furthermore, differentiating the equation (12.71) with respect to t and multiplying 
the result by d t Th give 

^\\d t T h (t, -,x 0 )\\ 2 L 2 + c 0 \\d t T h (t, -,x 0 )\\ 2 L 2 

< ^\\dtF h {t,-, xo)\\l 2 + {A h d t T h (t, ■,xo),d t r h (t, -,x 0 )) = 0, for t > 1, 
which further shows that 

\\dtTh(t, xo )\\ 2 L 2 < Ce c °^ 1 - ) •, ^o)||| 2 (( 1 / 4; i) ; i 2 (Qp < Ce Cot for t > 1. 

In a similar way one can derive ||c? tt F/ l (t, -,iCo)||i 2 < Ce~ Cot for t > 1. These 
inequalities together with (14.41) imply 

(4.23) || d t F(-, -,a;o)||Li(nx(i,oo)) + || •, xo)||Li(nx(i,oo)) < C, 

which together with (14.221) leads to (12.111) for the case h < ho := -RoFqj~ 2 /(16C*) 
with C * being given by (14.211) . 
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Secondly when h > ho, the decomposition in subsection 2.3 is not needed and 
the energy estimates of (12.6I) - (12.7[1 yield 

II d t Th\\L^Q T ) + \\dtT\\ L 2 iQT) < ||viyo)|| L 2 (n) + ||vr(o)|| L 2 (n) < ch ^^ /2 
||a«r h || LaWT) + ||ft t r|| L 2 (QT) < j|Vftr h (o)|U2 (n) + ||vd t r(o)|| L 2 (n) < ch~ 3 ~ N/2 , 
which imply 

(4.24) [ f (\d t F(t,x,Xo)\ + \tduF(t,x,Xo)\)dxdt <C. 

Jo Jn 

Since both 11 d t T h (t)\\ L 2 (n) +1| d tt T h (t)\\ L 2 ( n) an d ||9 t r(t)|| i 2 (n) + || l 9 tt r(t)|| i 2 (n) decay 
exponentially as t — > oo, it follows that (12.111) still holds when h > h 0 . 

The proof of Lemma 12.21 is completed. I 

5. Conclusion 

In this paper we have proved that the discrete elliptic operator — Ah gener¬ 
ates a bounded analytic semigroup and has the maximal L p regularity, uniformly 
with respect to h, in arbitrary convex polygons and polyhedra under the regularity 
assumption a tJ £ W 1 ' N+a . We have assumed the quasi-uniformity of the triangu¬ 
lation, and analysis of the problem under non-quasi-uniform triangulations remains 
open. As far as we know, only the analytic semigroup estimate m and its equiv¬ 
alent resolvent estimate were studied with an extra logarithmic factor for some 
special cases of non-quasi-uniform triangulations, see mm- The discrete maximal 
regularity estimates (11.101) and (11.131) have not been established with more general 
triangulations even in smooth settings. 

Appendix: The proof of Lemma 2.1 

Proof of \2.1\) : The inequality (12.11) is similar to Theorem 3.1.3.1 of [17], which 
was proved by using the local energy inequality of Lemma 3.1.3.2, and the lemma 
was proved under the assumption £ W' 1 , °°(f2), where Q is a convex domain. In 
the following, we show that this assumption can be relaxed to Oj 3 - £ VL 1 ,Ar+ “( 12 ) c —> 
C 7 (fl), where 7 = a/(N + a) 

Lemma A.l. If Q is convex and aij £ W l ' N+a (£i), then each point y £ Q has a 
neighborhood Bn(y) D such that 

(A.l) IMIipqn) < C\\Au\\ L 2^ + C||u||^i(n)- 

for all u £ H 2 D Hq such that the support of u is contained in Bn(y) D fl. The 
radius R depends only on the semi-norms |ay|wi,Ar+c(o) and \aij\ G1 fjiy 

Proof. Following the proof of Lemma 3.1.3.2 in [171 (see page 143, (3.1.3.4) and 
the equality above (3.1.3.5)), we have (using our notations) 

N 

IM| ff 2 (n) < C\\Au\\ L 2 ( Q) +C max \atj{y) - a lj (a;)|||u|| ff 2 (n) 

X G Vy 

i,j =1 

N 

H” C ^ ^ 
i,j =1 

< C'||Au|| i 2(f 2 ) + Ci?''||u||iy2(Q) 
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N 

+ C y~] ||Vay|| I ,iv+ a (n)||Vu|| I ,2(jv+ 0 .)/(jv-2+a)( n ). 
i,3 = 1 

When i? is small enough we have 

|M|ff 2 (n) < ^11^^11^2(0) + C'||Vu|| i 2(iV + c)/(JV-2 + Q )(0)- 

Since H 2 is compactly embedded into p{W 2 ( Ar+Q 0/( Ar-2+Q 0 which is again embedded 
into H 1 , there exists 8 a £ (0,1) such that 

||Vu|| i 2(N+c)/(jv— 2+c ,)(o) < e||M||ff 2 (Q) + C e ||w||Hi(o)> Ve G (0,1). 

Choosing e small enough, (IA.1I) follows from the last the last two inequalities This 
completes the proof of Lemma IA.1I I 

Then (12.11) can be proved by using Lemma fA. 1 1 and a perturbation procedure (as 
mentioned in the proof of [171 Theorem 3.1.3.1]). 


Proof of \2.2\) : Theorem 3.4 of [13] states that if C R w (TV > 3) is convex and the 
coefficients a,; ; are Holder continuous (so that (3.1)-(3.3) of [13] hold), the Green’s 
function of the elliptic operator A with the Dirichlet boundary condition satisfies 

(A.2) I VxG e {x,y)\ + \V y G e (x,y)\ < , 

where we have used the symmetry G e (x, y) = G e (y, x). Therefore, any Hq solution 
of the equation —V • («Vu) = / satisfies 


|Vu(s)| 


V x G e (x,y)f(y)dy 


< 


C 

\x — y\ N_1 


Il/llip(n) < C\\f\\ Lna) 

lp ' (n) 


for p > N. As pointed out in [T3] (page 227, the paragraph below Proposition 1), 
the inequality (1A.2I) for TV = 2 can be proved with some minor modifications on 
the proof of [T51 Theorem 3.3-3.4] since Theorem 3.3-3.4 of [T8J only requires ay- 
being Holder continuous coefficients. 


Proof of H2.3\) - j2.4\) : Since W 1 ’ N+a (fl) G(H), Theorem 1 of ]2I implies that 
the solution of the elliptic equation 

(A.3) { Au = 1 “ % 

y u = 0 on oil , 

with continuous coefficients a,j in a convex domain H C satisfies 

(A.4) llwllwc^o) < C'||/||w- 1 .<!(n)j Vl< 9 <oo, 

where we have noted that a continuous function is (5, R) vanishing, and a convex 
domain is (<5, a, i?)-quasiconvex [21] . Since the solution of (O) with f = 0 satisfies 
(integrating the equation against \u\ q ~ 2 u) 

IMIi'qn) < I|w°||l<i(o) ) 

it follows that the semigroup generated by the elliptic operator A is a contraction 
semigroup on L q (Q). By Theorem 1, Section 2, Chapter 3 of [43], the semigroup 
{E(t)}t >o generated by the elliptic operator A has an analytic continuation (ana- 
lyticity of the semigroup {£'(t)} t> o). Moreover, by the maximum principle we have 
u° > 0 =>• u > 0 (positivity of the semigroup {i?(t)}t>o) and then, by Corollary 
4.d of [50], the solution of m with u° = 0 has the maximal L p regularity (12.31) . 
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In other words, the map from / to Au given by the formula 

Au = f AE(t — s)/(-, s)ds 
Jo 

is bounded in L p ((0,T); L q ), for all 1 < p, q < oo. Since 

A 1/2 u = f A 1/2 E{t~s)f(-,s)ds = f AE(t- s)A- 1/2 f{-,s)ds, 
Jo Jo 


it follows that 


(A. 5) I|A 1/2 u|| L p ( (0,T);Z/i) < Cp,g\\A 1/2 /||lp((0,T);L«), VI < p,g < oo. 

It remains to prove the boundedness of the Riesz transform VA -1 / 2 : 

(A.6) IIWT 1 / 2 /!!^ < C q \\f\\ L <, V1 < g < oo. 

Then the last two inequalities imply 

||Vw||z,p((o,T);L9) = ||VA _1 / 2 ( J 4 1 / 2 u)|| i ,p((0 J T);I,a) 

< C'q||A 1 / 2 u|| I ,p(( 0)T ). i g) 

< C' Pig ||^l _1/2 /|| L p((0 i T);L?) 

< CpA/Ww-'-*, V1 < p, g < oo, 

where the last step of the inequality above is due to the following duality argument 
(At -1 / 2 j s self-adjoint): 

(A -1 / 2 /, g) = (f,A~ 1/2 g) < C\\f\\ w -i. 9 1| VA^gW^ < C\\f\\ w -x.,\\g\\ Lt >. 

It has been proved in [421 Theorem B] that the Riesz transform is bounded on 
L q {Vl) (i.e. the inequality (1A.6I) holds) if and only if the solution of the homogeneous 
equation 


(A.7) 


Au = 0 


satisfies the local estimate 

(A - 8) (wL , 

\' JnrB r (xo) J 




r(*o) 


|Vu| 2 diE 


1 

2 


for all xq € and 0 < r < ro, where ro and (7q > 2 are any given small positive 
constants such that f l fl B t7oro (x o) can be given by the intersection of B aoro (xo) 
with a Lipschitz graph. It remains to prove (1A.8[) . 

Let w be a smooth cut-off function which equals zero outside B 2r := B 2r {xo) 
and equals 1 on B r . Extend u to be zero on i? 2 r\^ and denote by U 2 r the average 
of u over B 2r • Then (IA.7I) implies 


N 

Y. di(aijdj(u(u - u 2r ))) 

i,j= 1 

N N 

(A.9) = ^2 di(a,ij(u — U 2 r)djUj) + o,ijdiwdj(u — u 2 r ) in 14, 

i,j =i »,j=i 


and the W 1,q estimate (IA.4I) implies 

||w(m - U2r)||iy 1 .®(fJ) < C\\(u - U 2 r)djW\\ Lq (n) + C\\diU}djU\\ W -l, q ^ 

< c\\(u - u 2r )djU}\\Li(n) + C\\diudju\\ L .(a) 
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= C\\(u - U2r)djU]\\ L q(B 2 r) + C\\diU)djU |U=(B 2r ) 

< Cv^WVuWls^), 

where s = qN/(q + N) < q satisfies L S (R) ^ W _1,9 (f2) and Vt rl ’ s (f2) ^ A«(R). 
The last inequality implies 

(A.10) IIVM|| L<! (onB r ) < C'r _1 ||Vu|| iS ((onB 2 r)- 

If s < 2 then one can derive 

l|Vii|| Li(nnB r ) < Cr N/q - N/2 \\Vu\\ L 2 (nnB2r) . 

by using one more Holder’s inequality on the right-hand side. Otherwise, one only 
needs a finite number of iterations of (1A.10I) to reduce s to be less than 2. This 
completes the proof of (1A.8[) . 

The proof of (12.3D - (I2.4D is complete. I 
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